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1.  INTRODUCTION 


This  paper  describes  a  study  which  used  elastic  wave  propaga¬ 
tion  theory  to  predict  and  analyze  ground  motions  produced  by  near  sur¬ 
face  airburst  explosions;  the  air-earth  environment  was  modeled  as 
three  homogeneous  elastic  layers  -  air,  soil  ana  rock  -  separated  by 
plane  parallel  boundaries  as  illustrated  in  Figure  1.  The  explosive 
source  is  located  on  the  axis  of  symmetry.  The  air  is  treated  as  an 
elastic  fluid,  while  the  soil  and  rock  are  treated  as  elastic  solids. 
Elastic  material  parameters  that  characterize  the  wave  propagation  are 
the  compress ional  wave  (P  wave)  speeds  ,  the  shear  wave  (S  wave) 

speeds  CSi  ,  and  the  densities  Ci  •’here  i  =  1,  2,  and  3  for  the 
air,  soil,  and  rock,  respectively. 


^ihe  exact  closed  form  integral  solutions  of  Cagniard  (1)  for 
the  reflection  and  refraction  of  spherical  waves  in  elastic  solids 
were  adapted  and  extended  to  model  the  ground  shock  propagation  in  a 
layered  earth.  In  this  formulation  the  particle  motion  is  obtained  as 
a  sum  of  components  propagated  along  rays  or  paths  (such  as  shown  in 
Figure  1)  associated  with  distinct  P  and  S  wave  arrivals.  Calcu¬ 
lations  using  the  Cagniard  procedure  were  used  previously  successfully 
by  the  author  to  predict  the  reflection  of  underwater  explosion  shock 
waves  from  the  ocean  bottom.  (See  References  (2)— (5> . )  The  theoreti¬ 
cal  analysis  and  computer  code  development  for  the  ground  shock  calcu¬ 
lations  were  extensions  of  the  bottom  reflection  study.  The  details 
of  the  theoretical  model  and  the  computer  code  used  in  the  elastic 
ground  motion  calculation  are  not  presented  but  will  be  given  in  a 
forthcoming  report  (6).  In  additio 
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Figure  1.  Model  for  airburst  explosions  over  layered  earth  media. 


provide  background  information  and  some  other  recent  applications  of 
the  Cagniard  elastic  propagation  theory. 


2.  MODELING  THE  AIRBLAST 


The  airblast  pressure  produced  by  an  explosion  attenuates 
with  radial  distance  R  from  the  source  more  like  L/R^  for  pressure 
levels  10  to  100  psi  than  the  spherical  elastic  1/R  decay  rate.  The 
blast  propagation  rate  decreases  with  range  instead  of  the  constant 
speed  Cpi  .  In  addition,  the  pressure-time  waveform  changes  in  shape 
with  range.  Thus,  the  airblast  cannot  be  directly  modeled  elastically 
but  must  be  approximated  by  either  (a)  linearizing  around  a  particular 
range  related  to  the  time  of  dominant  motion  or  (b)  by  simulating  the 
pressure  amplitude  and  arrival  time  by  a  distribution  of  sources.  The 
simpler  linearization  approach  was  used  in  this  study.  Several  proce¬ 
dures  were  investigated,  but  linearization  around  the  directly  trans¬ 
mitted  shear  wave  (path  A  of  Figure  1)  produced  the  best  agreement 


with  measured  waveforms  for  materials  ranging  from  weak  soils  to  hard 
rocks.  f 
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The  airblast  was  approximated  as  follows.  The  directly 
transmitted  shear  wave  path  A  of  Figure  1  was  determined  by  iterating 
for  the  range  rs  at  which  the  ray  enters  the  soil.  An  empirical 
formula  was  used  to  calculate  the  airblast  arrival  time  ta  for  an 
initial  value  of  rs  .  An  average  P  wave  speed  Cpi  was  computed 
from 


+  rs2)/ta 


Then  C  was  substituted  into  the  equation  for  Snell's  Law  of  acous¬ 
tics  toPobtain  a  new  estimate  for  rs  .  This  process  continued  until 
the  initial  and  final  values  of  rs  were  within  an  acceptable  toler¬ 
ance.  The  point  source  amplitude  and  pulse  shape  were  then  chosen 
using  empirical  formulae  so  that  the  airblast  pressure  was  matched  at 
the  point  (rg,  0). 

3.  COMPARISON  OF  CALCULATED  AND  MEASURED 
PARTICLE  VELOCITY  WAVEFORMS 


Calculations  were  performed  for  the  three  CENSE  (Coupling 
Efficiency  of  Near  Surface  Explosions)  explosive  field  test  series. 
(See  References  (10)  and  (11).)  These  tests  were  chosen  because  they 
provide  a  variety  of  site  characteristics  which  were  relatively  well 
controlled.  The  test  beds  were  either  effectively  homogeneous  or  had 
layering  suitable  for  the  two-laver  model.  In  addition,  each  of  the 
series  had  near  surface  airburst  explosions  for  which  particle  veloc¬ 
ity  or  acceleration  was  measured  in  the  upper  layer  for  a  variety  of 
ranges  from  the  explosions. 


CENSE  1  consisted  of  a  series  of  1000-lb  spheres  of  nitro- 
methane  detonated  over  a  massive  Kaventa  sandstone  formation.  These 
events  provide  data  for  cheeking  the  calculations  for  motion  in  a 
strong,  homogeneous  material  which  behaves  elastically  for  stress 
levels  of  hundreds  of  psi.  The  surface  rock  was  thick  enough  that 
three  layers  were  not  needed  for  the  computations.  Figure  2  compares 
the  theoretical  and  measured  vertical  and  radial  velocity  components. 
Vertical  velocity  is  positive  for  upward  motion  and  radial  velocity  is 
positive  for  outward  motion.  Note  that  the  experimental  and  theoreti¬ 
cal  curves  are  plotted  on  different  scales  and  that  the  calculation 
represents  only  part  of  the  measured  curve.  The  material  properties 
used  for  this  calculation  were  pj  ■  0.0012  gm/cm^  ,  Cp2  ■  9  ft/msec  , 
Cg2  ■  4  ft/msec  ,  and  c>2  “  2.4  gm/cm-*  .  Event  1,  shown  in  the  figure, 
was  detonated  with  its  charge  center  6  ft  above  the  rock.  Measurements 
were  made  with  velocity  gages  (having  a  nominal  600-Hz  frequency  re¬ 
sponse)  placed  2  ft  below  the  rock  surface.  The  airblast  peak 
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pressure  above  the  gages  was  measured  at  51  psi  at  the  48-ft  range. 
Other  calculations  not  shown  here  demonstrated  equally  good  agreement 
at  pressure  levels  from  10  to  120  psi. 

The  data  of  CENSE  2  provided  information  to  check  the  elas¬ 
tic  calculations  for  a  two-layered  clayey-silt  soil  site.  The  explo¬ 
sives  used  were  300-lb  spherical  TNT  charges.  In  the  calculations 
shown  in  Figures  3  and  4,  the  soil  was  modeled  in  two  layers:  a  sur¬ 
face  layer  20  ft  thick  with  Cp2  =  1.1  ft/msec  ,  Cs2  =0.6  ft/msec  , 
and  P2  =  1.7  gm/cra3  ,  and  a  lower  half-space  with  Cp3  =  1.6  ft/msec, 
Cs3  =  0.7  ft/msec  ,  and  p  =  1.75  gm/cm3  .  The  measurements  were  made 
with  velocity  gages  located  1.5  ft  below  the  surface.  Event  2  was 
detonated  with  charge  center  7.2  ft  above  the  soil  surface.  Excellent 
agreement  was  obtained  in  Figure  3  at  the  67-ft  (13-psi)  range.  There 
was  also  similarly  good  agreement  at  the  57-ft  (16-psi)  range.  Fig¬ 
ure  4  at  the  43-ft  (34-psi)  range  shows  slightly  poorer  agreement  but 
still  within  typical  scatter  of  field  measurements.  At  a  range  of 
32  ft  (60  psi)  the  linear  calculations  begin  to  fail  to  reproduce  the 
major  characteristics  of  the  measured  motion.  The  linear  theory  does 
not  predict  the  large  initial  downward  and  outward  displacements  seen 
in  the  experiments.  These  differences  are  probably  a  result  of  the 
nonlinear  material  properties  of  the  soil  becoming  important  and  a  re¬ 
sult  of  the  close-in  source  conditions  not  being  adequately  modeled 
by  the  localized  airblast  input  used  in  the  calculations. 

CENSE  3  provided  measurements  for  comparison  with  theory 
for  a  weak  soil  layer  over  a  hard  rock  site.  This  series  consisted 
of  seven  explosions  of  200  lb  (226  lb  TNT  equivalent)  of  nitromethane . 
The  test  bed  consisted  of  compacted  backfill  of  "alluvium"  soil  placed 
over  a  Kayenta  sandstone  deposit  similar  to  that  of  CENSE  1.  The 
thickness  of  the  soil  was  varied  from  0  to  6  ft.  Measurements  of  ver¬ 
tical  and  radial  acceleration  were  made  at  middepth  in  the  soil  layers 
and  in  the  rock.  Velocity  histories  were  obtained  by  integrating  the 
acceleration  records.  Events  2  and  4  were  surface  tangent  bursts, 
that  is,  the  explosive  charge  was  resting  on  the  soil  surface.  The 
soil  layer  thickness  in  Figure  5  was  6  ft.  In  Figure  6  the  thickness 
was  3  ft.  The  material  properties  used  in  the  calculations  were 
Cp2  =  0.9  ft/msec  ,  Cg2  =  0.3  ft/msec  ,  and  P2  ■  1.6  gm/cm^  for  the 
soil  and  Cp3  =  8  ft/msec  ,  Cs3  *  3  ft/msec  ,  and  P3  *  2.4  gm/cm3 
for  the  sandstone. 

At  the  56-ft  (12-psi)  range  of  Figure  5  the  calculations  are 
in  good  agreement  with  the  experimental  curves  up  to  a  time  of  about 
45  msec  if  the  high  frequency  spikes  are  neglected.  These  spikes  re¬ 
sult  from  using  an  airblast  pulse  with  zero  rise  time.  High  frequency 
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motion  of  this  type  is  filtered  out  of  the  measurements  because  of  the 
nonlinear  effects  of  the  soil  and  finite  frequency  response  of  the 
gages  and  recording  system.  At  the  32-ft  (52-psi)  range  of  Figure  6, 
agreement  is  poorer  but  the  initial  velocity  amplitudes  are  still 
close  to  the  measured  values.  At  a  range  of  24  ft  (110  psi)  the  calcu¬ 
lated  initial  peaks  are  nearly  a  factor  or  two  higher  than  the  experi¬ 
mental. 


Figure  7  is  presented  for  comparison  with  Figure  6  to  illus¬ 
trate  the  effect  of  increasing  the  soil  layer  thickness  at  the  CENSE  3 
site.  The  initial  portions  of  the  records  are  produced  by  the  directly 
transmitted  P  and  S  waves  and  are  not  dependent  on  the  soil  thick¬ 
ness.  The  later  motion  is  a  complicated  interaction  of  reflected 
waves  for  moderate  layer  thickness.  In  going  from  a  layer  thickness 
of  3  ft  as  in  Figure  5  to  the  6-ft  layer  in  Figure  6,  the  change  in 
frequency  of  the  motion  is  roughly  proportional  to  the  layer  thickness 
change,  but  the  waveforms  at  12-ft  thickness  do  not  follow  this 
pattern. 


From  the  few  waveforms  presented  here,  one  cannot  draw  gen¬ 
eral  conclusions  on  factors  affecting  the  period  and  amplitude  of  the 
low  frequency  motion.  A  detailed  parameter  study  and  analysis  will  be 
necessary  to  determine  how  the  motion  changes  in  going  from  very  thin 
to  very  thick  layers.  It  appears  that  simple  rules  of  thumb  based  on 
S  or  P  wave  layer  transit  times  will  be  valid  in  only  very  re¬ 
stricted  ranges  of  thickness  and  elastic  parameters. 

4.  CONCLUSIONS 

Prediction  of  velocity  waveforms  using  the  Cagniard  formula¬ 
tion  of  the  elastic  theory  and  the  localized  airblast  source  model  can 
be  expected  to  be  accurate  within  the  scatter  of  explosive  tests  mea¬ 
surements  for  times  up  to  about  one  cycle  of  the  low  frequency  motion 
and  for  airblast  overpressure  levels  at  the  gage  range  up  to  about  40 
psi  for  explosions  over  weak  soils  and  over  100  psi  for  strong  rocks. 
At  pressure  levels  in  the  range  of  40  to  100  psi  for  explosions  over 
soil,  the  elastic  theory  still  predicts  the  general  character  of  the 
motion  but  overestimates  the  peak  velocities  and  underestimates  the 
large  initial  downward  and  outward  displacements. 

Introduction  of  a  finite  rise  time  in  the  airblast  source 
pulse  is  desirable  to  eliminate  high  frequency  spikes  in  the  calcula¬ 
tions.  The  simple  localized  airblast  source  model  linearized  around 
the  directly  transmitted  shear  wave  may  be  a  major  contribution  to 
failure  of  the  calculations  at  higher  pressures  and  at  late  time  on 
waveforms. 
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The  linear  wave  propagation  model  can  produce  motion  wave¬ 
forms  for  homogeneous  sites  such  as  the  CENSE  1  sandstone  at  very 
small  computer  cost,  but  because  of  the  rapidly  increasing  effort  re¬ 
quired  to  calculate  the  late  time  portion  of  waveform  for  layered 
media,  routine  calculations  to  times  greater  than  about  three  shear 
wave  transit  times  of  the  layer  appear  to  be  more  expensive  than 
linear  finite  difference  methods  or  normal  modes  techniques.  Calcula¬ 
tions  with  the  Cagniard  theory  for  more  than  two  soil  layers  appear  to 
be  quite  expensive  except  for  early  time  motion  or  special  cases  where 
reflections  in  one  layer  can  be  neglected.  The  primary  applications 
for  computing  motion  waveforms  in  layered  media  appear  to  be  early 
time  motions  up  to  about  two  shear  transit  times  at  a  relatively  low 
computing  cost  and  minimal  effort  to  change  code  input  parameters. 

Since  the  theory  follows  rays,  the  composite  waveforms  can  be  dis¬ 
sected  to  study  the  contributions  of  individual  arrivals.  This  prop¬ 
erty  of  the  method  makes  it  ideal  for  studying  the  basic  characteris¬ 
tics  and  effects  of  the  controlling  parameters  of  wave  propagations  in 
layered  media. 
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